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ABSTRACT 

Impulsive solar energetic electrons are often observed in the interplanetary 
space near the Earth and have an attractive diagnostic potential for poorly un- 
derstood solar flare acceleration processes. We investigate the transport of solar 
flare energetic electrons in the heliospheric plasma to understand the role of trans- 
port to the observed onset and spectral properties of the impulsive solar electron 
events. The propagation of energetic electrons in solar wind plasma is simu- 
lated from the acceleration region at the Sun to the Earth, taking into account 
self-consistent generation and absorption of electrostatic electron plasma (Lang- 
muir) waves, effects of non-uniform plasma, collisions and Landau damping. The 
simulations suggest that the beam-driven plasma turbulence and the effects of 
solar wind density inhomogeneity play a crucial role and lead to the appearance 
of a) spectral break for a single power-law injected electron spectrum, with the 
spectrum flatter below the break, b) apparent early onset of low-energy electron 
injection, c) the apparent late maximum of low-energy electron injection. We 
show that the observed onsets, spectral flattening at low energies, and formation 
of a break energy at tens of keV is the direct manifestation of wave-particle in- 
teractions in non-uniform plasma of a single accelerated electron population with 
an initial power-law spectrum. 

Subject headings: Sun: flares - Sun: X-rays, gamma rays - Sun: activity -Sun: 
particle emission 



1. Introduction 

Solar flares are extremely efficient at accelerating electrons to non-thermal energies, 
which can subsequently be observed either by their emission at X-ray and radio wave- 
lengths or escaping along open magnetic field lines via direct electron measurements near 
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the Earth (see lAschwandenl 120021 ; iBrown fc Kontarll2005l; lBenal2008l. fo r a review). The first 
in-situ observations of energetic particles ( van Allen fc Krimigisl 1965 ) opened up the non- 



electromagnetic window of flare accelerated particle observations. Solar impulsive electron 
events detected in-situ genera lly display broken power-law energy distributions with lower 
energies having softer spectra (jLinlll985l ). These events also show near time-of-flight velocit y 
dispersion and a beamed pitch-angle distribution (e.g. lLinlll985l ; iKrucker et al.lll999l . 120071 ). 
From this evidence, it is often believed that such electrons propagate scatter-free from the 
Sun to the Earth (e.g. Wang et al.ll2006l ). The observed correlation between the spectral in- 



dices of energetic electron s at the Sun from X-ray data and the Earth from in-situ data (ILin 



19851 ; IKrucker et al.l 120071 ) is often viewed as an additional support for this model. At the 
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The standard model of Type III solar radio bursts (IGinzburg &; Zhelezniakovlll958l ) sug- 
gests that electron beams propagating in the ambient solar wind plasma from the Sun to 
the Eart h can excite L angmuir waves, which in turn generate escaping plasma radio emis- 
sion (see iMelrosd Il985l . for a review) . In a one-dimensional treatment assuming travel along 
magnetic field lines, as faster electrons outpace slower electrons a positive gradient is formed 
in velocity space. If this positive gradient gets large enough to start the generation of waves, 
electron energy is resonantly transferred to Langmuir waves in the backgroun d plasma. This 



transfer of energy reduces the g radient in velocity space forming a plateau (jVedenov et al. 



19621 ; iDrummond &; Pines! Il962j ). For a spatially limited electron beam cloud, the plasma 
waves generated at the front of the cloud are absorbed by the electron beam at the back 
allowing the electrons to travel through the corona with sma ll energy losses and with a ve- 



locity that decreases with time due to plasma inhomogeneity (jKontarll2001al ). Although this 
broad picture is often supported by observations, the detailed picture of electron transport 
and plasma radio emission is far from well-understood. This is largely due to electron beam 
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Solar impulsive electron events can span a broad range of energies, from a few keV 
to hundreds of keV (ILin et al.lll996l ). Their energy distribution forms a broken power- law 
spectrum with the break energy in deca-keV range. Despite often showing the near time- 
of-flight dispersion, lower energy ele ctrons appear to arrive sooner than expected from a 
scatter- free model (IWang et al.ll2006l ). Since the low energy electrons of a few keV should 
lose their energy collisionally in the low corona, these electrons are believed to be accelerated 



- 3- 



high in the corona (iLin et al.l Il996l ). Recent time-of-injection analysis (jWang et al.l 120061 ) 
assuming scatter-free propagation of solar energetic electrons suggests the existence of two 
electron populations, one low energy beam injected before the start of the type III burst and 
one high energy beam injected after the type III burst. 

In this Letter, we investigate the electron propagation from the Sun to the Earth taking 
into account the scattering of electrons by beam-driven plasma waves. We show, for the 
first time, that the generation and absorption of Langmuir waves by an electron beam in 
the non-uniform solar corona leads to the appearance of a break energy in the observed 
spectrum at the Earth and naturally explains the observed apparent early injection of low 
energy electrons. 



2. Modelling electron propagation in the heliosphere 



The transport of energeti c electrons in the heliospheric plasma is governed by a vari- 
ety of different processes (see iMelrosd Il990l . for a review). In this work we consider solar 
energetic electrons propagating along open magnetic field lines and assume their transport 
can be described one-dimensionally ignoring electromagnetic effects. Under this assump- 
tion, the evolution of the electron distribution function f(v,x,t) [electrons cm -4 s] and the 
spectral energy density of electron plasma waves W(v ,x,t) [ergs cm ~ 2 1 can be described 
self-consistently by the following kinetic equations (e.g. lKontarll2001al ) 
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dt ' dk dx dx dk 
where duj/dk = 3v^ e /v is the group velocity of Langmuir waves, k is the wavenumber of a 
Langmuir wave, j(v,x) = ^ £ ^ 2 f{ is the plasma wave growth rate, and 7 C and 7l are the 
collisional and Landau damping rates of waves respectively. The first term on the right hand 
side of both Equation (pQ) and (j2J) describes the resonant in t eracti on, a y = to, of electrons 



and L angmuir waves and was first derived by lVedenov et al.l (119621 ) and lDrummond fc Pines 
( 119621 ). W(v,x,t) is normalized to the wave energy density jWdk [ergs cm -3 ] and plays a 
similar role for plasma waves as the electron distribution function does for particles. 
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2.1. Initial electron beam distribution 

The initial distribution function is assumed to be a power-law with spectral index a in 
velocity space and has a finite spatial size d at initial time t = 0: 

f(v,x,t = 0) = g Q (v)exp (~^) ( 3 ) 

where 

So(«) = ™beam ( J , <* > 1 (4) 

is the initial electron distribution function normalized to nj, eam , the beam electron number 
density, v m i n is the low velocity cutoff, and a is the spectral index of the initial electron beam. 
The injected electron flux density differential in energy Fq(E, x, t = 0) [electrons cm -2 keV -1 
s _1 ] is also a power law Fq(E) ~ E~ 5 , where 5 = a/2. The initial spectral energy density 
of the Langmuir waves is assumed to be thermal W(v,x,t = 0) = /c j bT/(27t 2 A 2 :) ), where T 
is the background plasma temperature, ks is Boltzmann constant and A^e is the electron 
Debye length. 



2.2. Heliospheric plasma density 



The background p lasma is modelled u sing a heliospheric density model by iParkerl (119581 ) 
with normalization by iMann et al.l (119991 ) that can be found by numer ically integra ting the 
equations for a stationary spherical symmetric solution for solar wind (|Parkerlll958l ) 



r 2 n(r)v{r) = C = const 
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GM s /2v 2 , T is the electron plasma temperature, 



M s is the mass of the Sun. The constant appearing above is fixed by satellite measurements 
near the Earth's orbit (at r = 1 AU, n = 6.59 cm -3 ) and equates to 6.3 x 10 34 s -1 . 



3. Simulation of electron transport 



The system of kinetic equa tions ([U [2]) have been solved using finite difference methods 
as described in ( lKontarll2001al ) for a variety of initial beam parameters. The temperature 
of the heliospheric plasma was taken to be constant at T = 10 6 K and the plasma density 
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profile is given by the numerical solution of Equations (0 [6]). The low velocity cutoff was 
taken equal to approximately twice the thermal electron velocity v min = 1.2 x 10 9 cm/s. The 
initial spatial size of the electron cloud was taken as d = 5 x 10 9 cm so the injection time 
of electrons with velocity 5 x 10 9 cm /s is one se cond, which is a typical duration of type III 
bursts near the starting frequencies (jDulklll985l ). The initial beam density n^ eam was varied 
from 1 x 10~ 3 cm" 3 to 1 cm -3 to explain the observed fluxes near the Earth. The fluence 
spectra of such beam densities correspond to weak to medium solar energetic electron events 
observed near the Earth and are required for realistic simulation computational times. The 
initial height of the beam was 5 x 10 9 cm corresponding to the density 2.1 x 10 9 cm" 3 (local 
plasma frequency ~ 415 MHz). The initial beam spectral indices 5 = a/2 were between 2.5 
and 4.5, consistent with the observational values iKrucker et al.l (120091 ) . 



3.1. Electron spectra at 1AU 

Generally, it is seen from our simulations that as soon as the plasma wave growth time 
~ l/j(v, x) is less than the time scale of an electron cloud d/v for some energy E(v), the 
wave-particle interactions start to play an important role in the electron transport. Since this 
condition is energy dependent the electrons above certain energies are too dilute to generate 
plasma waves. The overall spectrum observed at the Earth becomes close to a broken 
power-law, where electrons below the break energy generate and absorb Langmuir waves, 
while electrons above the break energy are not affected by the wave-particle interactions. 



Traditionally in-situ measurements of energetic electrons (e.g. iLin et al.lll995l ) provide 
the flux density differential in energy F(E,x,t) = f(v,x,t)/m [electrons cm -2 s _1 keV" 1 ] 
and the fluences (flux integrated over the duration of an event) [electrons cm -2 keV" 1 ]. The 
injected electron fluence in our model J_ f(v,x,t)/vdx can be calculated from equations 
( l3f4"j) and is presented in Figure [H The corresponding energy spectral index of the injected 
electron fluence at the Sun is (a+l)/2 — 5 + 1/2. The resulting spectrum of solar energetic 
particles at the Earth is also presented in Figure [TJ As can be seen from Figure [TJ the 
spectrum of energetic particles above the break ~ 35 keV is identical to the spectrum of 
injected electrons so we can deduce these particles have indeed propagated scatter-free (in 
our model). The particles below the break energy do not propagate freely but generate 
electron plasma waves. The beam generated plasma waves dr ift in velocity- space toward 



lower phase velocities due to the solar wind density gradient (jKontarl l2001al ). This drift, 
caused by the decreasing ambient plasma density, takes waves out of resonance with the 
particles which generated them and so reduces the wave energy at a given point in phase 
space. At lower phase velocities these waves can be more easily absorbed by the thermal 
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plasma via Landau damping. Therefore, particles arriving later to this point in phase space 
are unable to restore the injected spectrum because they cannot absorb the same amount 
of energy from the waves. This results in a flatter energy spectrum of electrons below some 
energy where beam-plasma interactions are important (Figured]). 



3.2. Break energies and spectral indices 

We note the spectrum below the break is not exactly a power-law although resembles one 
closely, therefore to compare wit h observations, which have normally poor energy resolution 



AE ~ 0.1E (e.g. lLin et al.lll995l ). we fitted our simulated spectra with simple power-law fits. 
The spectral index below the break energy 8i ow is always smaller than the spectral index 
above the break energ y 8 high arid cor r elates (Figure [2J) with 6 high m a remarkably similar 



manner as observed by iKrucker et al.l (120091 ). The range of 8i ow appear in a rather narrow 
range between 2 and 2.5 for a wide range of injected spectral indices between 3 and 5 (Figure 
12]). The actual value of 8[ ow is also dependent on the background plasma density and will be 
different should the heliospheric density model change. 

The break energy range for all simulations is between 4 keV and 80 keV (Figure [2]), 
with the exact break energy being dependent on the initial spectral index of the beam, Shigh, 
and the initial density of the beam. As the initial density of the beam increases, the break 
of spectral index occurs at higher energies. Indeed, the larger number of injected electrons, 
the faster the generation of plasma waves proceeds and hence the stronger the interaction 
between electrons and plasma waves. This also explains the dependence of break energy to 
the injected spectral index, with lower spectral indices having a larger population of higher 
energy electrons and hence having higher break energies. The fluence at the break energy 
correlates to the break energy itself (Figure [2]) with higher break energies corresponding to 
lower fluence magnitudes. 



3.3. Electron time-of-flight and apparent injection time 



The particles arriving at the Earth (1.2 AU) show near time-of-flight dispersion (Fig- 
ure [3]), often observed by satellites in impulsive solar electron events. Assuming scatter-free 
propagation, i.e. without a ny interaction with p l asma, one can prod uce the apparent in- 
jection profile at the Sun, as IKrucker et al.l (119991 ) : IWang et al.l (120061 ) did for observations. 
Under scatter-free assumption, the time-of-flight for energetic electrons from the Sun to the 
spacecraft is simply = L/v(E), where L = 1.2 AU is the distance and v(E) is the speed of 
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Fig. 1. — The simulated spectrum (fluence [electrons cm -2 keV -1 ] ) of solar flare energetic 
particles at the Earth (crosses). The blue (red) line shows the power law fit to the spectrum 
below (above) the break energy 35 keV. The green line shows the initially injected fluence. 
Spectral index of injected electrons is 5hi g h = 4. The spectral index below the break is 
o low = 2.35. 
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Fig. 2. — Top: Spectral index above the break Shigh versus break energy for various electron 
densities and spectral indices; Middle: Spectral index below the break 5i ow energy versus 
spectral index above the break energy; Low: Fluence at the break energy versus break energy. 
Note that the simulations have been done for different particle densities and different spectral 
indices as given in the top panel. 
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electrons for various energies. The time of arrival is then simply tA = Unj + L/v(E), where 
t in j is the apparent injection time. These apparent injection profiles with background added 
have been calculated from the simulated fluxes at 1.2 AU (see Figure H]). If the electrons 
propagate scatter-free they would require 10 — 20 minutes earlier onset of injection t in j of low 
(3 — 12 keV) energy electron injection and a delayed maximum of the injection to explain 
observations. An identical simulation was ran with an electron beam not interacting with 
the background plasma (scatter free propagation) and the results are compared in Figure HI 
Comparing the calculated injection profiles over various energies one observes good agree- 
ment between free-streaming and full simulations at energies above the break energy 35 keV. 
The apparent injection profile at lower energies becomes wider than the free streaming case, 
and 3-25 keV electrons show apparent early injection of 1 - 20 minutes before the actual 
injection. The lower energy is, the earlier the apparent onset of the injection. The low 
energy electrons at the Earth are not only due to the injection at the Sun but because of 
the in-flight deceleration of faster particles. These electrons had initial energy higher than 
the detected energy at 1.2 AU, hence they have traveled a part of the distance faster than 
can be inferred from their detected energies. The relaxation of the electron distribution 
function towards a flatter shape in velocity space df(v,x,t)/dv ~ means at a specific 
spatial location, some electrons have energies too low to have arrived by free propagation 
alone. In addition, the maximum of apparent injection profile at energies below ~ 35 keV 
appears later. Therefore, the similar injection profile obtained (jWang et al.l 120061 ) should 
be interpreted as the direct evidence of electron plasma wave scattering in the heliosphere 
and not the indication of a separate acceleration mechanism. This early injection time is a 
direct result of low-energy electron driven turbulence in non-uniform plasma, which affects 
the propagation of electrons. As evident from Figure HI the onset of electron injection is 
also instrument background dependent - the higher/lower background levels would lead to 
later/earlier injection times for low energy electrons. 



4. Discussion and Conclusions 

The generation and reabsorption of electron plasma waves by the electron beam in a non- 
uniform plasma plays an important role in the electron transport and should be taken into 
account when in-situ electron measurements are analyzed and interpreted. The simulations 
presented successfully reproduce the spectral and temporal characteristics of observed solar 
energetic electron events. The scattering of the solar energetic particles by the beam-driven 
electrostatic plasma waves leads to the appearance of a broken power-law in energy spectrum, 
and the apparent early injection of low energy electrons in a few keV range at the Sun. We 
emphasize here that it is the combined effect of plasma wave generation and non-uniform 
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Fig. 3. — Simulated elect ron flux density time profiles of energetic electrons for WIND/3DP 
energies ( iLin et al.lll995l ). Electron flux density [electrons cm" 2 s" 1 keV" 1 ] as a function of 
time at 1.2 AU for 10 energy channels. The time t — corresponds to the injection time 
at the Sun. The sawtooth structure appearing in low energy channels is an artifact of finite 
binning in the velocity space. The initial beam parameters are the same as in Figure [H 
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Fig. 4. — The apparent injection profiles of electrons at the Sun assuming free streaming of 
all electrons. The plotted fluxes are normalized to the maxima of the scatter free case. The 
injection profiles for scatter-free propagation (without generation and absorption of waves) 
is overplotted with dashed lines. The initial beam parameters are the same as in Figure [3J 
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density inhomogeneity that leads to the appearance of a broken power-law. 

We have shown that the low energy electrons (below the break) are originally injected 
with higher energies but have lost some their energy to plasma waves in the background 
plasma and are therefore detected earlier than their energy at the spacecraft suggest. The 
apparent early start of low energy electron injection a ppears due to pro pagation effects 
and does not require a secondary beam population (c.f. IWang et al.ll2006l ). Moreover, the 
simulations naturally explain the linear correlation between low energy an d high energy 



spect ral indices often observed in the spectra of solar impulsive electron events (IKrucker et al. 
2009h . 



The characteristic time of beam-plasma interaction via electron plasma waves is in- 
versely proportional to the density of the energetic electrons. If the beam is dilute, electrons 
do not generate plasma turbulence and the spectrum of such electrons could be free from 
propagation effects (in our model). Such events are likely to be seen only at low energies (e.g. 
Krucker et al.ll2009l ). If the beam is dense enough to excite Langmuir waves, the initially 
injected power-law spectrum will be detected as a broken power-law. The break energy is 
dependent on a number of parameters: spectral index of injected solar electrons, the density 
of the energetic electrons, and the heliospheric density model. Therefore, the correlation 
between break energy and fluence at the break energy should be taken with care. Our 
simulations also suggest that in order to deduce meaningful time-of-flight parameters (e.g. 
injection times), the energetic electrons should be always analyzed in the channels above the 
break energy and not above some fixed energy as sometimes done in the literature. 

While our simulations can successfully explain and reproduce a number of observed 
properties and highlight the governing role of wave-particle interactions at keV- tens of 
keV energies, the heliospher ic plasma - electron beam syste m is more complex. First, the 
different density models (e.g. lNewkirklll967l ; ISaito et al.lll977l ) could change the break energy 
values, although will not alter the overall conclusions. Second, the heliospheric plasma 
has density perturbation s on various scales which can affect the evolution of plasma waves 



(IKontar fc Pecseli 



Kontar 2001b 



Li et al 



2002), lead to a spiky structure of the Langmuir waves (IMelrose 11990 



20061 ). often observed in the interplanetary space (jLinl Il985l ). and 



hence affect the electron distribution. In addition, the higher energy solar impulsive electron 
events demonstrate injection with ~ 10 min delay after the onset of type III bursts (e.g. 
Krucker et al.lll999l ; lHaggerty fc Roelolll2002l ). Therefore, additional simulations and in-situ 
measurements are needed to understand this complex non-linear system. 



EPK acknowledges financial support of STFC rolling grant and STFC/PPARC Ad- 
vanced Fellowship. Financial support by the European Commission through the SOLAIRE 



13 



Network (MTRN-CT-2006-035484) is gratefully acknowledged. 

REFERENCES 

Aschwanden, M. J. 2002, Space Science Reviews, 101, 1 

Benz, A. O. 2008, Living Reviews in Solar Physics, 5, 1 

Brown, J. C, & Kontar, E. P. 2005, Advances in Space Research, 35, 1675 

Cane, H. V. 2003, ApJ, 598, 1403 

Drummond, W. E., & Pines, D. 1962, Nucl. Fusion SuppL, 3, 1049 
Dulk, G. A. 1985, ARA&A, 23, 169 

Ergun, R. E., Larson, D., Lin, R. P., McFadden, J. P., Carlson, C. W., Anderson, K. A., 
Muschietti, L., McCarthy, M., Parks, G. K., Reme, H., Bosqued, J. M., D'Uston, C, 
Sanderson, T. R., Wenzel, K. P., Kaiser, M., Lepping, R. P., Bale, S. D., Kellogg, P., 
& Bougeret, J.-L. 1998, ApJ, 503, 435 

Gaelzer, R., Ziebell, L. F., Vifias, A. F., Yoon, P. H., & Ryu, C.-M. 2008, ApJ, 677, 676 

Ginzburg, V. L., & Zhelezniakov, V. V. 1958, Soviet Astronomy, 2, 653 

Gosling, J. T., Skoug, R. M., & McComas, D. J. 2003, Geophys. Res. Lett., 30, 130000 

Grognard, R. J.-M. 1982, Sol. Phys., 81, 173 

Haggerty, D. K., & Roelof, E. C. 2002, ApJ, 579, 841 

Kontar, E. P. 2001a, Sol. Phys., 202, 131 

-. 2001b, A&A, 375, 629 

Kontar, E. P., & Pecseli, H. L. 2002, Phys. Rev. E, 65, 066408 

Krucker, S., Kontar, E. P., Christe, S., & Lin, R. P. 2007, ApJ, 663, L109 

Krucker, S., Larson, D. E., Lin, R. P., & Thompson, B. J. 1999, ApJ, 519, 864 

Krucker, S., Oakley, P. H., & Lin, R. P. 2009, ApJ, 691, 806 

Ledenev, V. C, Zverev, E. A., & Starygin, A. P. 2004, Sol. Phys., 222, 299 



-14- 



Li, B., Robinson, P. A., & Cairns, I. H. 2006, Physical Review Letters, 96, 145005 
Lin, R. P. 1985, Sol. Phys., 100, 537 

Lin, R. P., Anderson, K. A., Ashford, S., Carlson, C, Curtis, D., Ergun, R., Larson, D., 
McFadden, J., McCarthy, M., Parks, G. K., Reme, H., Bosqued, J. M., Coutelier, J., 
Cotin, F., D'Uston, C, Wenzel, K.-P., Sanderson, T. R., Henrion, J., Ronnet, J. C, 
& Paschmann, G. 1995, Space Science Reviews, 71, 125 

Lin, R. P., Larson, D., McFadden, J., Carlson, C. W., Ergun, R. E., Anderson, K. A., 
Ashford, S., McCarthy, M., Parks, G. K., Reme, H., Bosqued, J. M., d'Uston, C, 
Sanderson, T. R., & Wenzel, K. P. 1996, Geophys. Res. Lett., 23, 1211 

Lin, R. P., Potter, D. W., Gurnett, D. A., & Scarf, F. L. 1981, ApJ, 251, 364 

Mann, C, Jansen, F., MacDowall, R. J., Kaiser, M. L., & Stone, R. G. 1999, A&A, 348, 614 

Melnik, V. N., Lapshin, V., & Kontar, E. 1999, Sol. Phys., 184, 353 

Melrose, D. B. 1985, Plasma emission mechanisms (Solar Radiophysics: Studies of Emission 
from the Sun at Metre Wavelengths), 177-210 

-. 1990, Sol. Phys., 130, 3 

Newkirk, G. J. 1967, ARA&A, 5, 213 

Parker, E. N. 1958, ApJ, 128, 664 

Saito, K., Poland, A. I., & Munro, R. H. 1977, Sol. Phys., 55, 121 

van Allen, J. A., & Krimigis, S. M. 1965, J. Geophys. Res., 70, 5737 

Vedenov, A. A., Lelikhov, E. P., & Sagdeev, R. Z. 1962, Nucl. Fusion Suppl., 2, 465 

Wang, L., Lin, R. P., Krucker, S., & Gosling, J. T. 2006, Geophys. Res. Lett., 33, 3106 

Zaitsev, V. V., Mityakov, N. A., & Rapoport, V. O. 1972, Sol. Phys., 24, 444 



This preprint was prepared with the AAS IATgX macros v5.2. 



